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ABSTRACT 

We analyze in detail the hydrodynamics and afterglow emission of an ultrarelativistic blast wave when 
it expands in a density-jump medium. Such a medium is likely to appear in the vicinity of gamma-ray 
bursts (GRBs) associated with massive stars. The interaction of the blast wave with this medium is 
described through a reverse shock and a forward shock. We show that the reverse shock is initially 
relativistic if the factor of a density jump (a) is much larger than 21, and Newtonian if 1 < a <C 21. 
We also calculate light curves of the afterglow emission during the interaction if the reverse shock is 
relativistic, and find that the optical flux density initially decays abruptly, then rises rapidly, and finally 
fades based on a power-law, which could be followed by an abrupt decay when the reverse shock has 
just crossed the originally swept-up matter. Therefore, one property of an afterglow occurring in a 
large-density-jump medium is an abrupt drop followed by a bump in the light curve and thus provides 
a probe of circumburst environments. In addition, this property could not only account for the optical 
afterglows of GRB 970508 and GRB 000301C but also explain the X-ray afterglow of GRB 981226. 

Subject headings: gamma-rays: bursts — relativity — shock waves 



1. INTRODUCTION 

The density distribution of circumburst environments is 
one of the most important issues in the theories of gamma- 
ray bursts (GRBs) (for review articles see Piran 1999; van 
Paradijs, Kouveliotou & Wijers 2000; Meszaros 2001). On 
one hand, it is directly pertinent to the progenitors of 
GRBs. Two currently popular models for the progenitors 
are the mergers of compact stars (neutron stars or black 
holes) and the explosions of massive stars. It has been 
argued that GRBs produced by the former model occur 
in a uniform interstellar medium (ISM) with density of 
~ 1 cni^'^ and GRBs in the latter model occur in pre-burst 
winds (Chevalier & Li 1999) and/or giant molecular clouds 
(Galama & Wijers 2001; Reichart & Price 2001). Thus, 
an environmental signature may provide a clue about the 
GRB progenitors. On the other hand, the environmental 
properties can directly influence the decay rates of after- 
glows. For example, afterglows arising from the interac- 
tion with pre-burst winds should decay more rapidly than 
afterglows do in a low-density medium (e.g., ISM) (Dai 
& Lu 1998a; Meszaros, Rees & Wijers 1998; Panaitescu, 
Meszaros & Rees 1998; Chevalier & Li 1999, 2000). Fur- 
thermore, ultrarelativistic fireballs (or jets) in a uniform 
dense medium (e.g, galactic-like giant molecular clouds) 
must evolve to the non-relativistic regime within a few 
days after the bursts, leading to a rapid decay of the af- 
terglows (Dai & Lu 1999, 2000; Wang, Dai & Lu 2000). It 
is thus natural that an afterglow signature can probe the 
ambient matter as well as the progenitors. 

In the previous afterglow shock models, the environ- 
ments of GRBs are usually assumed to be continuous 
media (e.g., ISM and wind). Actually, there are possi- 
bly jumps (or bumps) in the density profile of the am- 
bient media of GRBs associated with massive stars. Such 
jumps may be produced by several astrophysical processes, 
e.g., the deceleration of winds in their external medium 



(Ramirez-Ruiz et al. 2001; R. Wijers 2001, private com- 
munication) or the interaction of fast and slow winds (Luo 
& McCray 1991; Vikram & Balick 1998) or pre-burst su- 
pernova ejecta (Vietri & Stella 1998). In this Letter, we 
perform a careful analysis for the afterglow emission when 
a post-burst relativistic blast wave interacts with such a 
density-jump medium. We first analyze the hydrodynam- 
ics of the interaction in §2 and then discuss the afterglow 
signature in detail in §3. Our findings are summarized and 
discussed in S4. 



2. HYDRODYNAMICS 

Let's envision that some central energy source produces 
an ultrarelativistic fireball. After the internal shock emis- 
sion (viz., a GRB), the fireball will start to sweep up its 
ambient medium, leading to an ultrarelativistic blast wave. 
The initial hydrodynamics of this interaction has been an- 
alyzed in detail by Sari & Piran (1995), and studied numer- 
ically by Kobayashi, Piran & Sari (1999) and Kobayashi 
& Sari (2000). Now we consider this ultrarelativistic blast 
wave which first expands in an ISM or a stellar wind and 
then hits an outer high-density region. We assume that the 
medium of interest has a simple density profile: n — AR~'' 
for R < Rq and n = rti = constant for R > Rq. Here 
A — uq X 1 cm^^ if the inner medium is an ISM (s = 0), 
and A = 3 X lO'^^A, cm^^ if the inner medium is a wind 
(s = 2). Such a density profile seems to be able to recon- 
cile the contradiction that the multiwavelength afterglow 
fits in the jet model by Panaitescu & Kumar (2001) indi- 
cate an ambient medium density significantly lower than 
that expected in star-forming regions. Ramirez-Ruiz et al. 
(2001) showed that the interaction of winds from massive 
stars with the ambient medium can indeed lead to a large 
density jump at a few lO^'' cm. 

Before the blast wave hits the high-density medium, 
its Lorentz factor decays based on the Blandford-McKee's 



2 



(1976) solution as 7 = 8.2£;5^^no ^^^t'^/^ in the inner- 

ISM case and 7 = S-SE^^A'^^S''^/'^ in the inner-wind 
case, where £^53 is the energy of the blast wave in units 
of 10^^ ergs and t is the observer's time in units of 1 day 
(neglecting the redshift effect). Subsequently, the inter- 
action of the blast wave with the high-density medium is 
described through two shocks: a reverse shock that prop- 
agates into the hot shell (viz., the R < Rq medium swept 
up by the blast wave), and a forward shock that propa- 
gates into the high-density mcdimn. Therefore, there are 
four regions separated by the two shocks in this system: 
(1) unshocked high-density medium, (2) forward-shocked 
high-density mcdiiun, (3) reverse-shocked hot shell, and 
(4) unshocked hot shell. We denote rij, Cj and pi as the 
baryon number density, energy density and pressure of re- 
gion "i" in its own rest frame respectively; ji and /3, are 
the Lorentz factor and dimensionless velocity of region "i" 
measured in the local medium's rest frame respectively; 
and 7ij and /3y are the relative Lorentz factor and dimen- 
sionless velocity of region "i" measured in the rest frame 
of region "j" respectively. If 7^ ^ 1 and jj ^ 1, then 
jij ~ {ji/jj +7i/7i)/2- We further assume the equations 
of state for regions 2, 3 and 4 to be relativistic and region 
1 to be cold. Thus, the equations describing the jump con- 
ditions for the forward and reverse shocks become (Bland- 
ford & McKee 1976; Sari & Piran 1995; Kumar & Piran 
2000; Zhang & Meszaros 2001a) 
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= 72-1, 



n2 



= 472 + 3, 



(1) 



In this case, the reverse shock converts only a small frac- 
tion ([^[ <C 1) of the kinetic energy into thermal energy 
because 274 ~ y/f. Thus, the forward shock expands al- 
most at the velocity of the previous blast wave. 

In the next section, we will discuss light curves of the 
afterglow emission if the reverse shock is initially relativis- 
tic. For this purpose, we need to know how the thermody- 
namic quantities and Lorentz factor of each region evolve 
with radius R at two different stages: 

When the reverse shock crosses region 4, this region al- 
ways expands adiabatically at a constant Lorentz factor 
of 74, and thus we have rii oc R~^ and 64 oc nl^^ oc R~^. 
As a result, we obtain / oc R~'^. For regions 2 and 3, 
we have 72 = 73 oc /^/^ oc R~^, 62 = 63 oc 7I oc i?^^, 
712 oc 72 oc R~^, and ns oc ^4 (63/64) oc R~^. At this 
stage, 734 oc R, which is different from the initial hydro- 
dynamics of a relativistic blast wave, 734 oc R^/^, derived 
from Sari & Piran (1995) in the thin shell case. This is 
due to a hot region 4 as compared to Sari & Piran (1995). 

After the reverse shock crosses region 4, the profile of 
the shocked high-density medium begins to approach the 
Blandford-McKee solution as long as the shocked matter 
has a relativistic equation of state, as shown numerically 
at the initially hydrodynamic stage of an afterglow by 
Kobayashi et al. (1999) and Kobayashi & Sari (2000). Us- 
ing this solution, the Lorentz factor and the energy density 
of a given fluid element for region 3 decay as 73 oc 
and 63 oc i?-26/3 (Sari & Piran 1999a, 1999b). In addi- 
tion, the Lorentz factor and the energy density of region 2 
evolve as 72 oc and 62 oc R~^. 



73^ 



(l + 363/64)(3 + 63/64) 
1663/64 

_ (63/64)(l -I- 363/64) 
3-1-63/64 



(2) 



(3) 



where rUp is the proton mass. 

Regions 2 and 3 should keep the pressure equilib- 
rium and velocity equality along the contact discontinuity, 
which yield 72 = 73 and 62 = 63. Under these conditions, 
the solution of equations (l)-(3) depends only on two pa- 
rameters: 74 and / = 64/(nimpC^). The solution has two 
limits which correspond to the cases that the reverse shock 
is relativistic or Newtonian. If 63 ^ 64, then the reverse 
shock is initially relativistic: 



^1/2^1/4 31/4^^2 

72 = 73 = 31/4 ^ 734 = -JJTJ4- > 1> 



(4) 



which requires a = n\/no » 64/3 ~ 21, where no is the 
baryon number density of the inner medium at R = Rq and 
the energy density of region 4 at this radius is assumed to 
equal A'y'lnQmpC^. From equation (4), 73 <C 74, showing 
that most of the initial kinetic energy of region 4 is con- 
verted into thermal energy by the shocks. On the other 
hand, for 1 < a <C 21, the reverse shock is Newtonian: 



734 - 1 ^ 



274/v7-i y_i 



274/\/7+2/\/3y 2 



or 



72 = 73 ~ 74(1 - 



(5) 



(6) 



3. LIGHT CURVES OF THE EMISSION 

We now consider synchrotron radiation from all the re- 
gions at two different stages. The electron energy dis- 
tribution just behind the shock is usually a power-law: 
due/dje oc 7~^ for 7e > 7m- Here we discuss only 
the case of p > 2. Dai & Cheng (2001) have discussed 
light curves of the emission from a relativistic shock in 
an ISM or a wind for 1 < p < 2. Assuming that Cg 
and eB are constant fractions of the internal energy den- 
sity going into the electrons and the magnetic field re- 
spectively, we have the electron minimum Lorentz factor, 
7m,i = [{P - 2)/(p - l)]eeei/(niTOeC^), and the magnetic 
field, Bi = {SneBeiY^^, for region "t", where irie is the 
electron mass. 

According to Sari, Piran & Narayan (1998), the spec- 
trum consists of four power-law parts with three break 
frequencies: the self-absorption frequency, the typical syn- 
chrotron frequency i'm,i — lilm i^-^i / i'^'^^ec) , and the 
cooling frequency = l^Trem^c/ {a^j.Bf^it^), where ctt 
is the Thomson scattering cross section. In this Letter 
we neglect the self-absorption because it does not affect 
the optical radiation which we are interested in. In or- 
der to calculate the flux density at a fixed frequency, one 
still needs to derive the peak flux density. The observed 
peak flux density is given by F^^^i = Ne^i'yiP^^^i/{A7rDj^), 
where Ng^i is the electron number of region "i" at radius R, 
Pvm.i = n^eC^o'TBi/{3e) is the radiated power per electron 
per unit frequency in the frame comoving with the shocked 
matter, and Dl is the source's luminosity distance to the 
observer. So we in fact need to calculate A^e.i- First, it 
is easy to obtain the total electron number of region 2 
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by N(.,2 — (47r/3)ni(_R'' — R^). Second, the time interval 
for the reverse shock to spend in crossing a length inter- 
val dx'i in the rest frame of region 4 is dt'^ = dx'^/{(334C), 
which is related to a time interval {dt';^) in the rest frame 
of region 3 by dt'^ = 734rft'j, and furthermore dt'^ is re- 
lated to a time interval (dim) in the local medium's rest 
frame by dt'^ = = dR/{'ys/3sc), so we calculate the 

electron number of region 3: Ne.s ~ J^^'^ A-KR^fudx'^ = 
/^^47ri?2n4[734/334/(73/33)]di? (X - i?2)^ ^here (i^^ « 
f3z K, 1, and thus the electron number of region 4 is 
iVe.4 = 47rAi?Q^''/(3 — .s) — A^e,3, when the reverse shock 
crosses region 4. Letting iVe,3 equal the initial electron 
number of region 4 (viz., when 7Ve,4 = 0), we can calculate 
the radius at which the reverse shock has just crossed this 

region: i?A = -Ro(l + 5.83 x lO'^ny^n^y^) for s = and 
Ra = Ro{l + 1-25 X 10-'^E~^^^^ A^n^l^^to^^^) for s = 2, 
where ni^s = ni/10^cm~^ and to is the observer's time (in 
1 day) at i? = Rq. 

Although the spectrum does not depend on the hydro- 
dynamics of the shocked matter, the light curve at a fixed 
frequency is determined by the temporal evolution of ^^,1, 
fc,< and These quantities depend on how ji, nt, ej 

and Ne^i scale as a function of R as well as t either for 
-Ro ^ -R < Ra or for R > Ra- However, wc note that for 
the typical values of the involved parameters (see below), 
Ra/Ro — 1 "C 1, showing that 7i, rii and are almost un- 
changed and thus we consider only the temporal evolution 
of i^cA (fx t~^) and Ne^i for Ro < R < Ra- 

One crucial effect in calculating the light curve is that 
the photons which are radiated from different regions at 
the same time measured in the local medium's rest frame 
will be detected at different observer times (Zhang & 
Mcszaros 2001a). The understanding of this effect is that 
for Rq < R < i?A the Lorentz factor of region 4 is much 
larger than that of regions 2 and 3 so that for a same time 
interval in the local medium's rest frame, dR/c, the emis- 
sion from region 4 reaches the observer in a time interval 
of ^ dR/ {2c^1), while the emission from regions 3 and 2 
reaches the observer in a time interval of ~ dR/{2c^^). 
After considering this effect and the scaling law of 7^ with 
R, wc obtain the observer time of the radiation from re- 
gion 4 at i? = Ra- tA,4 = to{l -1-1.16 x lO^^ny^n^g^^) for 
s = and tA,4 = io(l + 2.50 x 10"2£;-^i/2^^^-i/2^-i/2^ 
for s = 2, and the corresponding observer time of the ra- 
diation from regions 3 and 2: ^a.s = 1.33fo for s = and 
tA,3 = 2.0*0 for s = 2, where we have used Ra/Rq — '^ 1. 
It is thus seen that tA.i/ta — 1 <C f-A.^/ta — 1, implying that 
the observed radiation from regions 3 and 2 is indeed de- 
layed as compared with that from region 4. 

Fig. 1 presents two R-band {vr ~ 4.4 x 10^^ Hz) 
light curves about the effect of an ultrarelativistic blast 
wave interacting with a density-jump medium on the af- 
terglow. The outer-medium density assumed here, rii ~ 
10"^ cm^'', is the one of typical galactic-like giant molecu- 
lar clouds. It can be seen from this figure that at t > to 
the flux density (F^„ ) initially drops abruptly, then rises 
rapidly, and finally declines based on a power-law fol- 
lowed by an abrupt decay at t = tA,3- This result is 
easily understood: an initially abrupt decay of the emis- 
sion is due to the spectral cutoff frequency fcut,4 < t^R 



{solid line) or the rapid decrease of the electron num- 
ber m region 4, Ne,4 cx [AtAA/ta - (t/ta + if] {dashed 
line), during the adiabatic expansion for to < t < tA,4- 
Note that VcutA = i^c,o{Rl Ro)~'^ results from evolution 
of the fast-cooling electrons of region 4, where Vc.o = 

2.8 X IQ^'^e^^-yEll^'^n^^tQ^^'^ Hz is the cooling frequency of 
the originally swept-up electrons sX t = to- In the period 

01 to < t < tAAi the emission flux densities from regions 

2 and 3 are low both because the two shocks have swept 
up only a small number of the electrons and because the 
radiation from these regions reaches the observer at a later 
time than the radiation from region 4 does. As the number 
of the electrons swept up by the two shocks increases, the 
flux density increases rapidly as ^Vj^ oc {t/to — l){t/to)~^ , 
where the first factor arises from Fy^^i oc {t/to — 1) and 
the second factor from Vc.i oc {t/to)~'^ for regions 2 and 3. 
However, since i^cs < t^m,3 < J^r for the parameters shown 
in the figure, all the electrons in region 3 are in the fast 
cooling regime. As a result, the R-band flux density of the 
radiation from this region disappears at t > tA,3 and thus 
only the radiation from region 2 could be detected. There- 
fore, an abrupt decay of the flux density could appear at 
t = tA,3- We further define the factor of emission brighten- 
ing (TZ) as the ratio of the observed density fluxes with and 
without a density jump at t = tA,3- In Fig. 1, 7?. ~ 5.1 
{solid line) and TZ ^ 10 {dashed line). Fig. 2 c;xhibits 
the light curves of the afterglow emission when the inner 
medium at i? < i?o is a stellar wind. In this flgure, TZ ~ 5.3 
{solid line) and 7?. ^ 11 {dashed line)- Recently Ramirez- 
Ruiz et al. (2001) also estimated TZ ~ a^P+i)/!**^-**)! ~ 4.5 
and 21 for s = and s = 2 (where p = 2.5 and a = 10^) re- 
spectively. However, we find that TZ depends on the shock 
parameters (e.g, Cg and es) for a fixed density jump. One 
reason for the difference between our and Ramirez-Ruiz 
et al.'s (2001) results is that for to < t < tA.3 the optical 
emission arises partially from the reverse shock considered 
here. Another reason is that the electrons which produce 
the optical emission may be in different radiation regimes 
before and when the reverse shock crosses region 4. 

4. DISCUSSION AND CONCLUSIONS 

Wc have analyzed the hydrodynamics and afterglow 
emission of an ultrarelativistic blast wave when it inter- 
acts with a density-jump medium. This interaction is de- 
scribed through two shocks: a reverse shock and a forward 
shock. We have shown that the reverse shock is initially 
relativistic if a 21, and Newtonian if 1 < a <C 21. We 
have also investigated in detail light curves of the afterglow 
emission during the interaction if the reverse shock is rel- 
ativistic, and found that the R-band flux density initially 
decays abruptly, then rises rapidly, and finally fades based 
on a power-law, which could be followed by an abrupt de- 
cay when the reverse shock has just crossed the originally 
swept-up matter. 

Our analysis is basc;d on several simplifications: First, 
we have considered only one density jump, but our dis- 
cussion should be in principle applied to a more realistic 
case in which there could be a few jumps (or bumps) in 
the density profile. Second, since the Y parameter for syn- 
chrotron self-Compton (SSC) is not far larger than unity 
in our model, the effect of SSC on the optical emission is 
insignificant. This effect was recently discussed by several 
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authors (e.g., Waxman 1997; Wei & Lu 1998; Panaitescu & 
Kumar 2000; Sari & Esin 2001; Zhang & Meszaros 2001b). 
Third, we have assumed a spherical relativistic blast wave, 
but an actual blast wave may be a jet, whose edge effect 
and sideways expansion can lead to a steepening of the 
hght curve (Meszaros & Rees 1999; Rhoads 1999; Sari, Pi- 
ran & Halpern 1999; Dai & Cheng 2001). Fourth, we have 
neglected the emission off the line of sight. Such emission 
could significantly alleviate the initial abrupt drop of the 
flux density. 

The humps have been observed to appear in the light 
curves of several optical afterglows (e.g., GRB 970508 and 
GRB 000301C). In our model, these humps are under- 
stood to be due to the contribution of the radiation from 
regions 2 and 3 when the reverse shock crosses region 4. 
A few other interpretations have been proposed, e.g., re- 
freshed shocks due to Poynting-flux-dominated or kinetic- 
energy-dominated injection (Dai & Lu 1998b; Panaitescu, 
Meszaros k Rees 1998; Zhang & Meszaros 2001a, 2001c; 
Chang et al. 2001), and microlensing events (Garnavich, 
Loeb & Stanek 2000). However, there is not any abrupt 
drop in the afterglow light curves of the latter models. 
This property could be used to distinguish between the 
present model and the other interpretations. We will carry 
out detailed fits to the optical afterglows of GRB 970508 



and GRB 000301C based on the present model. In addi- 
tion, this model can account well for the rise-decline fea- 
ture of the X-ray afterglow light curve of GRB 981226 
(Frontera et al. 2000). 

Finally, we have noted that when the reverse shock 
crosses region 4, neutrinos with energies of TeV-PeV are 
possibly produced by ir^ created in interactions between 
accelerated protons and synchrotron photons from accel- 
erated electrons in regions 2 and 3. Furthermore, such 
neutrino emission is delayed about tA,3 after the GRB. 
This result is different from the prompt neutrino emis- 
sion discussed by many authors, e.g, Waxman & Bah- 
call (1997, 2000), BahcaU & Meszaros (2000), Meszaros 
& Rees (2000), Dai & Lu (2001), and Meszaros & Wax- 
man (2001). One expects that these delayed neutrinos, if 
detected, could provide further diagnostics about circum- 
burst density-jump environments. 
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Fig. 1. — R-band light curves of the afterglow emission when an ultrarelativistic blast wave interacts with a density-jump medium. The 
blast wave expands within an ISM (s = 0) until it reaches a high-density medium at Rq = 5 x 10^^ cm (or to = 0.26 days). The model 
parameters are taken; E53 = 1, no = 1, ni,3 = 1, te = 0.1, p = 2.5 and = 2 x 10^* cm. The solid and dashed lines correspond to ts = 0.1 
and 0.01 respectively. 
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